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We present the results of a systematic study of the shape of the pion distribution in coordinate 
space at freeze-out in Au+Au collisions at RHIC using two-pion Hanbury Brown- Twiss (HBT) 
interferometry. Oscillations of the extracted HBT radii vs. emission angle indicate sources elongated 
perpendicular to the reaction plane. The results indicate that the pressure and expansion time of 
the collision system are not sufficient to completely quench its initial shape. 

PACS numbers: 25.75.Gz, 25.75.Ld 



Relativistic heavy ion collisions arc believed to reach 
sufBciently high energy densities and temperatures for 
the possible formation of a quark-gluon plasma (QGP) 
[11. Hanbury Brown-Twiss (HBT) interferometry of 



two particle Bosc-Einstcin correlations directly accesses 
the space-time structure of the emitting source formed 
in these collisions, providing crucial probes of the sys- 
tem dynamics. At the Relativistic Heavy Ion Collider 
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(RHIC), identical-pion HBT studies in Au+Au collisions 
at y/I^N = 130 GeV 0, 0] yielded an apparent source 
size consistent with measurements at lower energies, in 
contrast topredictions of larger sources based on QGP 
formation |5|. In addition, hydrodynamical models, suc- 
cessful at RHIC in describing transverse momentum spec- 
tra and elliptic flow , have failed to reproduce the small 
HBT radii 0. This so-called "HBT puzzle" might 
arise because the system's lifetime is shorter than pre- 
dicted by models. 

In non-central collisions, azimuthally-sensitive HBT 
measurements performed relative to the reaction plane 
provide a measure of the source shape at freeze-out 0, 
llll H^ . In such collisions, the almond-shaped col- 
lision geometry generates greater transverse pressure 
gradients in the reaction plane than perpendicular to 
it. This leads to stronger in-plane expansion (elliptic 
flow) |l[ll[ll[l3 which diminishes the initial out-oj- 
plane spatial anisotropy. Therefore the freeze-out source 
shape should be sensitive to the evolution of the pres- 
sure gradients and the system lifetime; a long-lived sys- 
tem would be less out-of-plane extended and perhaps in- 
plane extended. Hydrodynamic calculations |l6l | predict 
a strong sensitivity of the HBT parameters to the early 
conditions in the collision system and show that, while 
the system may still be out-of-plane extended after hy- 
drodynamic evolution, a subsequent rescattering phase 
[l7j tends to make the final source in-plane. Knowledge 
of the freeze-out source shape might discriminate among 
scenarios of the system's evolution. 

In this Letter, we present results of a systematic study 
of azimuthally-sensitive HBT in Au-|-Au collisions at 
y/sNN = 200 GeV. These results allow for first studies of 
the relationship between the initial and final eccentrici- 
ties of the system. 

The measurements were made using the STAR detec- 
tor at RHIC. Particle trajectories and momenta were 
reconstructed using a Time Projection Chamber (TPC) 
with full azimuthal coverage, located inside a 0.5 Tesla 
solenoidal magnet. Au-|-Au events with primary ver- 
tices < 25 cm longitudinally of the TPC center were 
placed into centrality classes following Ref. A high- 
multiplicity triggered dataset of 500k events was used 
for the most central bin (0-5% total cross section), and a 
minimum bias dataset of 1.6 million events was used for 
all other centrality classes (5-10%, 10-20%, 20-30% and 
30-80%). The second-order event plane angle ^'2 
for each event was determined from the weighted sum 
of primary charged-particle transverse momenta |2ll |. 
Within a resolution which we determine from the random 
subevent method , ^2 ~ ^rp (true reaction plane an- 
gle) or ^"2 ~ + tt; i.e. the direction of the impact 
parameter vector is determined up to a sign |2(1 . . 

Pion candidates, selected according to their specific 
energy loss (dE/dx) in the TPC in the rapidity range 
\y\ < 0.5, were required to pass within 3 cm of the pri- 



mary vertex and contain > 15 (out of 45) TPC space 
points in the reconstructed trajectory. Pion pairs were 
subjected to two requirements. To account for recon- 
structing a single particle trajectory as two tracks, a 
topological cut is applied in which a minimum fraction of 
TPC pad layers must show distinct hits for both tracks. 
To reduce the effect of merging two particle trajectories 
into a single reconstructed track, an additional topolog- 
ical cut requires that the number of merged TPC hits 
falls below a maximum fraction. The latter cut leads to 
a systematic error that depends on the event multiplicity 
and the transverse momentum of the tracks . 

Pairs of like-sign pions were placed into bins of $' = 
'/'pair ~ ^2 J whcrc (/)pair is the azimuthal angle of the pair 
momentum [k = i(pi + P2)]. Because we use the 2"^^- 
order reaction plane, <&' is only defined in the range (0, tt). 
For each bin, a three-dimensional correlation function is 
constructed in the Pratt-Bertsch "out-side-long" decom- 
position of the relative pair momentum q. The nu- 
merator of the correlation function contains pairs of pions 
from the same event, and the denominator contains pairs 
of pions from different events which have similar primary 
vertex position, reaction plane orientation, multiplicity, 
and magnetic field orientation, pairs and 7r+ pairs 
were mixed separately due to charge-dependent accep- 
tances but are combined to increase statistics; separate 
7r+ and tt^ analyses showed no significant differences. 

Finite reaction plane resolution and finite width of the 
$' bins reduce the measured oscillation amplitudes of 
HBT radii vs. A model-independent correction pro- 
cedure 131, applied to each q-bin in the numerator and 
denominator of each correlation function, accounts for 
these effects and increases the amplitudes of the HBT 

The increase is roughly 



radu vs. <I> ($ 



inversely proportional to the measured |20l l21j reaction 
plane resolution, i.e. the amplitudes increase ~ 10-30%. 
All data were corrected using this procedure. Also, auto- 
correlation contributions to $ were tested by selecting 
distinct sets of particles for event plane determination 
and HBT analysis, with no observed effect. 

In addition, correlations due to final-state Coulomb 
repulsion must be accounted for, in order to isolate the 
Bose-Einstein correlations of interest. Traditionally this 
was accomplished by applying correction weights (deter- 
mined by calculating the Coulomb correlation function 
A'(q) for a spherical Gaussian source 0) to all pairs 
in the denominator. Recently, the CERES collabora- 
tion I23 noted that this approach over-corrects for the 
Coulomb effect and advocated an improved procedure 
[2^ which applies the Coulomb weight only to the frac- 
tion of pairs that participate in the Bose-Einstein correla- 
tion. We adopt this approach, and fit each experimental 
correlation function to the form: 

C{<l,^)=N- [(l-A)-l + A./^(q)(l + G(q,$))], (1) 

where the (1 — A) and A terms account for the non- 
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FIG. 1: (color online) Squared HBT radii using Eq.Qrelative 
to the reaction plane angle for three centrality classes. The 
solid lines show allowed l24l fits to the individual oscillations. 



participating and participating fractions of pairs, respec- 
tively, A'^ is a normalization parameter, and (^(q, <&) is 
the Gaussian correlation model 



G(q,$) 



(2) 



i?? arc the squared HBT radii, where the Z,s,o subscripts 
indicate the long (parallel to beam), side (perpendicular 
to beam and total pair momentum) and out (perpendic- 
ular to qi and qs) decomposition of q with an additional 
cross terixL Fitting with Eq.^ caused i?o to increase 
10-20% compared to Coulomb-correcting all pairs, while 
i?s and i?/, respectively, are consistent within errors. 

Figuren]shows the squared HBT radii, obtained using 
Eq. n as a function of <f> for three centrality classes. All 
pairs with pair transverse momentum 0.15 < < 0.6 
GeV/c are included, and each centrality is divided into 
12 $ bins of 15° width. The data point at $ = tt is the 
reflected $ = value, and solid lines indicate Fourier 
expansions of the allowed oscillations : 



(i?2(fcT,$)cos(n$)) (/i = o,s,/) 



(i?^(fcT, $) sin(ri,$)) (fi^^os) 



(3) 



As expected , the 0*'^-order Fourier coefRcient (FC) in- 
dicates larger apparent source sizes for more central col- 
lisions. We verified that the 0*'^-order FC corresponds to 
the HBT radii from an azimuthally-integrated analysis. 

Strong 2'^'^-order oscillations are observed for R^, i?^ 
and Rig , and the signs of the oscillations are quahtatively 
self-consistent [lol |24|. though the amplitude for most- 
central events is small. Similar oscillations were observed 
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FIG. 2: (color online) Squared HBT radii relative to reac- 
tion plane angle for four kr (GeV/c) bins, 20-30% centrality 
events. The solid lines show allowed | 2^ fits to the individual 
oscillations. 



in a statistics-limited analysis of minimum-bias Au+Au 
collisions at yijvw = 130 GeV [i^. These oscillations 
correspond to a pion source spatially extended perpen- 
dicular to the reaction plane, as discussed below. The 
next terms (4*'^-order) in the Fourier expansions (Eq. O 
are consistent with zero within statistical errors. 

The feT-dependence of the oscillations of the HBT radii 
may contain important information on the initial condi- 
tions and equation of state of the system ■ Figure [21 
shows the $ dependence of HBT radii for mid-central 
(20-30%) events for four kr bins. Due to the additional 
division of pairs in kx, only four bins in $ are used. 
The 0*'^-order FC increases with decreasing kx, which 
was observed for azimuthally-integrated HBT analyses 
at y^SNN = 130 GeV Q and attributed to pion emission 
from an expanding source. Strong out-of-plane oscilla- 
tions are observed for all transverse radii in each fcy bin. 

The full results are summarized in Figure |31 which 
shows the centrality dependence of the Fourier coeffi- 
cients for three ranges of fcy. The number of partici- 
pants for each centrality was determined using a sim- 
ple nuclear overlap model Systematic variations of 
the HBT radii arise due to their sensitivity to the anti- 
merging cut threshold and uncertainty associated with 
the Coulomb procedure ■ The total variation is largest 
for Rl Q (~ 10%). The systematic variation on the rela- 
tive amplitudes plotted in the right panels of Fig. |2| are 
negligible compared to statistical errors. Also, all corre- 
lation functions composing Fig. 13 are corrected for mo- 
mentum resolution following our prescription in Ref. 01 ■ 
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FIG. 3; (color online) Fourier coefBcients of azimuthal oscil- 
lations of HBT radii vs. number of participating nucleons, for 
three (GeV/c) bins. Left panels: means (0"'-order FC) of 
oscillations; right panels: relative amplitudes (see text for de- 
tails). Larger participant numbers correspond to more central 
collisions. 



As in Figs. □ and [3 the O^'-order FCs (left panels) 
correspond to the squared HBT radii that would be ob- 
tained in a standard analysis. R^o^ -^s o ^^^'^ ^lo 
all observed to decrease for more peripheral collisions. 
Ro/ Rs, found in theoretical calculations to be sensitive to 
the emission duration of the system Q , is observed to be 
Rofi/Rsfi = 1.15±0.01 (1.06±0.01) for the lowest (high- 
est) bin for 0-5% most-central events. These values 
are consistent with that reported at ^Jsnn = 130 GeV Q 
when the increase in Ro due to the improved Coulomb 
correction (Eq.QJ is accounted for. Ro/Rs is still smaller 
than the predictions from hydrodynamical models, indi- 
cating the "HBT puzzle" persists at ^/snn = 200 GeV. 
Dynamical effects on the homogeneity region affect 
as well as i?^ ^ [3, E3- 

The relative ampli- 
tudes of the oscillations offer a more robust measure of 
the spatial anisotropy and arc less sensitive to dynamical 
effects [s^ . FigureElshows (right panels) the relative am- 
plitudes vs. number of participants for three fc^ ranges, 
using the ratios Ra,2/^l,o (o^ = o,s,os) and Rf2/Rfo- 
The relative amplitudes for all three transverse radii de- 
crease in magnitude with increasing number of partici- 
pants, and their weak fc^ dependence agrees qualitatively 
with hydrodynamic calculations |l6l |. 

To extract the shape of the pion source at freeze-out, 
a model-dependent approach is required. In the pres- 




initial 



FIG. 4: Source eccentricity obtained with azimuthally- 
sensitive HBT (ennai) vs. initial eccentricity from a Glauber 
model (einitiai). The most peripheral collisions correspond 
to the largest eccentricity. The dashed line indicates einitiai = 
Eflnai- Uncertainties on the precise nature of space-momentum 
correlations lead to 30% systematic errors on Sflnai H^j. 



ence of collective flow the HBT radii correspond to re- 
gions of homogeneity [sJl and do not reflect the entire 
source. The "blast- wave" parametrization [il Isol Is^ Is^ 
of freeze-out, which incorporates both spatial and dy- 
namical anisotropics, has been used to describe various 
observables at ^/snn = 130 GeV [sO, 13- ^ recent 
blast-wave analysis yOl showed that the relative oscil- 
lation amplitudes {e.g. shown in Fig. 3) are most sen- 
sitive to the spatial anisotropy. The source eccentricity 
(e = {Rl - Rl)/{Rl + RD) can be related to the relative 
amp litude of the HBT oscillations by efinai ~ 2i?^ 
where R,, (Ry) is the radius of the elliptical 
source in-plane (out-of-plane). 

The eccentricity of the initial almond-shaped overlap 
region was calculated from a Glauber model ,19] using 



the r.m.s. values for Ry and Rx. Figure 0] shows the 
relation between the initial and final eccentricities ob- 
tained by averaging the three fc^ bins in Fig. |21 The 
initial and final eccentricities exhibit a monotonic rela- 
tionship, with more peripheral collisions showing a larger 
final anisotropy. Within this model-dependent picture, 
the source at freeze-out still retains some of its initial 
shape, indicating that the outward pressure and/or ex- 
pansion time was not sufficient to quench the initial spa- 
tial anisotropy. The large elliptic flow and small HBT 
radii observed at RHIC energies might favor a large 
pressure build-up in a short-lived system. Also, out- 
of-plane freeze-out shapes tend to disfavor a long-lived 
hadronic rescattering phase following hydrodynamic ex- 
pansion [T^ . 

In conclusion, we have performed an analysis of two- 
pion HBT interferometry relative to the reaction plane 
in Au+Au collisions at ^Jsnn = 200 GeV. The rela- 
tive amplitudes of the HBT radius oscillation is largest 
for peripheral collisions, indicating larger out-of-plane 
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anisotropy in the pion source at freeze-out, for colli- 
sions with larger initial spatial anisotropy. No strong kT 
dependence of the relative oscillation amplitudes is ob- 
served. The out-of-plane freeze-out shape of the source 
indicates that the build-up of pressure and the evolu- 
tion time of the expanding system are not sufficient to 
quench the initial geometry of the collision. This infor- 
mation, taken together with the size of the source and 
anisotropics in momentum space, places significant con- 
straints on future theoretical efforts to describe the na- 
ture and timescalc of the collision's evolution. 
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